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ABSTRACT 


The  purpose  of  this  paper  is  to  describe  and  discuss  results  obtained 
from  a  new  stress  routine,  which  is  implemented  in  the  Projectile  Design 
Analysis  System  (PRODAS).  This  system  is  regularly  used  by  the  Aerodynamics 
Branch  (FXA),  Air  Force  Armament  Laboratory  (AFATL),  Eglin  Air  Force  Base, 
Florida,  to  design  projectiles  and  to  predict  the  aerodynamic  behavior  and 
performance  of  projectiles  and  rockets  prior  to  testing  in  the  Aeroballistic 
Research  Facility  (ARF).  Due  to  the  increasing  variety  of  aerodynamic 
configurations  that  are  being  tested,  a  stress  routine  was  desired  which 
provides  the  model/sabot  designer  with  a  good  estimation  of  the  projectile 
stress  during  the  launch  phase. 
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SECTION  I 


INTRODUCTION 

The  free  flight  ballistic  range  has  been  and  still  is  an  important  tool 
in  the  testing  and  development  of  military  weapon  systems  and  ammunition  as 
well  as  ballistic  research.  The  design  of  the  models  used  in  ballistic 
research  is  a  critical  element  in  this  process.  As  a  result,  a  computer 
program  entitled  "The  Projectile  Design  and  Analysis  System"  (PRODAS) 
(Fig.1)  was  developed  to  predict  the  mass  properties  along  with  the 
aerodynamic  parameters  and  the  associated  free  flight  behavior  of  various 
munitions.  Another  purpose  of  this  program  is  to  provide  the  design 
engineer  with  information  about  the  expected  acceleration  loads  during  the 
launch  cycle.  For  a  complete  description  of  the  PRODAS  program  refer  to 
Ref.  1. 

PRODAS  is  basically  a  design  tool  which,  as  described  herein,  uses  the 
interactive  and  display  capabilities  of  the  Eglin  Air  Force  Base  graphic 
computer  system.  The  previously  existing  stress  routine  in  PRODAS  only 
performed  a  stress  analysis  for  spin  stabilized  projectile  configurations. 
The  new  stress  routine  discussed  herein  computes  the  projectile  compressive 
loads  in  the  outer  shell  at  numerous  locations  along  the  projectile  length. 
For  certain  projectile  geometries  containing  a  cylindrical  body,  the 
dynamic  material  properties  are  taken  into  consideration  and  the  maximum 
allowable  dynamic  stresses  are  compared  to  the  computed  compressive 
stresses  in  the  body.  The  basis  for  these  dynamic  properties  are  empirical 
results,  taken  from  several  sources  (Ref.  2-4). 

The  purpose  of  this  paper  is  to  discuss  this  new  stress  routine  and  to 
present  some  typical  results.  However,  it  should  be  noted  that  the 
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addition  of  this  new  stress  routine  does  not  represent  the  final  solution  to 
providing  the  engineer  with  design  information.  It  is  expected  that  this 
routine  will  be  further  improved  and  a  finite  element  routine  is  already  in 
development  which  will  also  be  included  in  the  near  future.  This  finite 
element  routine  may  be  the  subject  of  a  future  paper. 
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SECTION  II 


LAUNCH  INDUCED  STRESS 


1.  General  Case 

The  maximum  force  (F),  acting  on  a  projectile-base  during  the 
acceleration  phase  is  equal  to  the  maximum  base,  pressure  (P) 
delivered  by  the  propellant  multiplied  by  the  area  (A)  of  the  bore 
(Ref.  5,  6) 

F  =  P  A  (1) 

The  setback  force,  caused  by  the  acceleration  (a)  of  the  projectile  is 

F  =  a  W  /  g  (2) 

p 

where  g  is  the  acceleration  of  the  gravity  in  ft/sec  ,  and  W  is 
the  total  weight  of  the  projectile  in  pounds.  Combining  equations 
( 1 )  and  ( 2)  leads  to 

a  =  P  A  g  /  W  (3) 

The  inertia  of  the  mass  of  the  parts  of  the  projectile  ahead  of  a 
transverse  section  will  lead  to  a  compression  force  (FQ)  in  that 
particular  cross-section,  assuming  the  projectile  is  acting  as  a 
rigid  body. 

Fc  =  W'  a  /  g  (4) 

W'  is  the  weight  of  all  projectile  parts  forward  of  the  transverse 
section.  The  compressive  stress  is  then  defined  as: 

cr c  =  Fc 7  Ai  (5) 

Where  A^  is  the  cross  sectional  area  of  the  load  carrying  transverse 
section.  This  approach  is  applicable  to  projectile-models  with  a  rigid 
body  or  with  thick  shell  walls. 
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In  the  case  of  a  rifled  gun,  a  tangential  force  (Ffc)  also  exists 
which  is  caused  by  the  angular  acceleration  (a')  imparted  by  the 
rotating  band  on  the  shell.  This  angular  acceleration  is  a  function 
of  the  rifling  twist  (n,  in  calibers  per  turn),  the  linear 
acceleration  and  the  projectile  diameter  (d,  in  inches). 

a '=  24  JT  a  /  (n  d)  (6) 

The  torque  applied  to  the  projectile  is 

T  =  a ’  ( I/g )  (7) 

where  I  is  the  polar  moment  of  inertia  of  the  projectile  (lb. in.2) 
and  T  has  the  units  of  lb. in.  The  tangential  force  can  be  written  as 

Ft  =  T  /  (d/2)  (8) 

Combining  equations  (3),  (6),  (7)  and  (8),  we  can  obtain 


F 


t 


48  7T  I  P  A 

n  d2  W 


(9) 


Equation  9  shows  that  Ffc  is  directly  proportional  to  the  propellant 
pressure  acting  on  the  base  of  the  model  and  therefore  Ffc  will  be 
a  maximum,  when  the  base  pressure  is  a  maximum.  The  previously  existing 


PRODAS  program  contained  a  simplified  analysis  which  considered  the  shear 


stress  caused  by  the  tangential  force  (Ft)  applied  by  the  rotating  band. 

It  should  be  noted  that  there  are  other  forces  which  can  contribute  to 


the  stress  levels  experienced  during  launch.  For  example,  any  projectile 
with  internal  cavities  containing  a  filler  material  (i.e.  a  high  explosive 
HEI  round)  can  have  longitudinal,  tangential,  and  radial  stresses  resulting 
from  the  rotation,  setback,  or  movement  of  filler  material  or  any  other 
internal  components.  As  mentioned  previously  the  purpose  of  the  present 
work  was  to  incorporate  an  additional  routine  in  PRODAS  where  the 
compressive  stresses  acting  on  a  rigid  or  semi-rigid  projectile  are  caused 
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by  the  setback  forces  encountered  during  launch  (see  equation  5).  This  is 
applicable  for  saboted  projectiles  fired  from  a  smooth  bore  gun  where  the 
shear  stresses  resulting  from  rotation  are  negligible. 

In  order  for  the  designer  to  determine  whether  or  not  the  existing 
compressive  stresses  are  high  enough  to  possibly  cause  failure,  they  must  be 
compared  with  the  maximum  allowable  stress.  Since  the  existing  compressive 
loads  are  applied  in  a  dynamic  manner  and  exi3t  only  for  a  short  period  of 
time  (i.e.  Microseconds)  the  maximum  allowable  dynamic  load  (Q)  can  be 
significantly  higher  than  the  maximum  allowable  static  load  (see  Ref.  6). 


This  maximum  allowable  dynamic  load  Q  (lb.)  can  be  calculated 
by  the  following  secant  formula: 


Where  m  is  normally  set  equal  to  1.7  and  L  =  length  of  column 
(in.),  r  =  least  radius  of  gyration  of  column  section  (in.),  E  = 
modulus  of  elasticity  (psi),  A  =  section  area  of  column  (in.  ),  and  Cy  is 
the  static  yield  stress  (psi). 

Since  this  equation  is  nonlinear  in  Q,  it  can  only  be  solved  by  trial 
and  error  or  by  the  use  of  prepared  charts  (see  Refs.  7,8  and  the  attached 
appendix) . 

Under  certain  circumstances  the  maximum  load  a  body  will  sustain  is  not 
given  by  the  strength  of  the  material,  but  by  the  stiffness  of  the  body. 
This  behavior  is  known  as  "elastic  stability"  and  arises  when  the  load 
produces  a  bending  or  a  twisting  moment  that  is  proportional  to  the 
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corresponding  deformation.  An  example  of  this  is  the  Euler  column,  which 
is  a  straight  column,  axially  loaded.  It  remains  straight  and  suffers  only 
axial  compressive  deformation  under  small  loads.  If  while  thus  loaded  it 
is  slightly  deflected  by  a  transverse  force,  it  will  straighten  after 
removal  of  this  force.  But  there  is  some  axial  "critical  load"  that  will 
hold  the  column  in  the  deflected  position,  and  since  both  the  bending 
moment  due  to  the  load  and  the  resisting  moment  due  to  the  stress  are 
directly  proportional  to  the  deflection,  the  load  required  to  hold  the 
column  in  the  defected  state  is  independent  of  the  amount  of  the 
deflection.  Any  increase  in  the  "critical  load",  leads  immediately  to  a 
collapse  of  the  column. 

A  very  thorough  discussion  of  the  general  problem,  with  detailed 
solutions  of  many  cases  are  given  in  Ref.  6  and  7,  from  which  many  of  the 
formulas  presented  in  the  Appendix  were  taken. 

2.  Special  Case 

A  special  model  case,  which  is  representative  of  many  of  the  subscale 
models  tested  in  the  Aeroballistic  Research  Facility  (ARF),  was  defined  as 
follows:  The  model  has  a  cylindrical  body,  with  two  concentric  holes, 
drilled  from  the  base  of  the  projectile  towards  the  tip.  The  model  may 
consist  of  two  different  materials  where  the  nose  section  and  the  body 
section  is  joined  with  a  threaded  stud.  This  threaded  stud  can  be  either 
part  of  the  nose  section  or  the  body  section.  The  nose  section  may  also 
consist  of  various  elements  such  as  an  ogive  and  conical  elements  capped 
with  a  hemispherical  nose  tip  (see  Fig.  2). 
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For  the  above  defined  projectile,  the  previously  discussed  stress 
analysis  is  performed  and  then  compared  with  the  maximum  allowable  dynamic 
stresses  as  calculated  at  both  the  base  and  joint.  If  the  projectile  does 
not  fit  the  special  case  as  defined  above,  only  the  general  stress  analysis 
will  be  computed  and  the  design  engineer  will  be  left  to  his  own  means  in 
determining  whether  or  not  the  calculated  stresses  are  critical. 
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SECTION  III 


RESULTS 

When  running  the  new  stress  routine  in  PRODAS,  the  program  will 
automatically  determine  weather  or  not  the  conditions  for  the  "special 
case"  projectile  exist.  A  projectile  design  will  be  treated  as  follows: 

The  standard  stress  analysis  corresponding  to  the  previously  discussed 
method  will  be  computed  for  that  projectile.  The  stress  will  be  calculated 
at  200  transverse  sections,  beginning  at  the  projectile  tip  and  ending  at 
the  projectile  base.  The  longitudinal  distance  from  one  transverse  section 
to  the  next  is  equal.  The  information  about  the  acceleration  is  taken  from 
the  PRODAS  interior  ballistic  routine,  and/or  can  be  chosen  by  the 
designer.  Results  appear  in  the  form  of  tables,  as  shown  in  Table  1,  and 
plotted  versus  the  projectile  length,  as  shown  in  Figs.  3  and  4.  These 
results  provide  the  design  engineer  with  the  opportunity  to  redesign  that 
specific  model,  for  instance  in  the  joint  area  to  avoid  inappropriate  stress 
concentrations . 

In  addition  to  the  above  mentioned  stress  analysis,  the  maximum 
allowable  dynamic  stress  will  also  be  calculated  if  the  conditions  for  the 
specially  defined  projectile  exist.  In  order  for  this  to  be  accomplished 
it  is  necessary  for  the  designer  to  choose  the  materials  used.  This 
selection  is  made  from  the  table  as  shown  in  Table  2.  Depending  on  what 
materials  are  selected,  subtables  will  appear  on  the  screen  for  the 
designer  to  specify  certain  material  properties  (i.e.  the  maximum  yeild 
point)  of  the  selected  material. 

"Enter  the  yield  strength  of  the  material 
(cylindrical  part)  in  10^  psi.  To  keep 
the  default  value  of  68  (hit  'return')" 
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The  computed  maximum  allowable  dynamic  stresses  are  then  displayed  for  both 
the  base  and  joint  cross  sections  of  the  specially  defined  projectile. 


Projectile  Base: 

Stress:  9791  psi 

Dynamic  Allowed  Stress:  39*217  p3i 
Safety  Margin:  4.00 

Joint  Area: 

Stress:  9106  psi 

Dynamic  Allowed  Stress:  49,848  psi 
Safety  Margin:  5.47 
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SECTION  IV 


CONCLUSIONS 

A  Fortran  V  subroutine  has  been  included  in  the  Projectile  Design  and 
Analysis  System  (PRODAS)  in  order  to  analyze  the  compressive  stress  along 
a  projectile  body  during  launch.  Also,  the  maximum  allowable  dynamic 
stresses  are  computed  for  a  specially  defined  projectile.  It  is  believed 
that  this  new  stress  routine  will  be  of  great  assistance  to  the  design 
engineers  of  the  Aeroballistic  Research  Facility  and  will  significantly 
reduce  the  risk  of  launch  failures  due  to  inadequately  designed  models.  It 
is  expected  that  this  routine  will  be  further  improved  in  the  future  (i.e. 
by  adding  a  sabot  analysis)  and  that  more  advanced  routines  (i.e.  finite 
element)  will  also  be  incorporated. 
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Figure  1:  PRODAS  STRESS  MENU 
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?ib,le  1 Stress  analysis  results  for  the  'standard  projectile' 
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4348.3 

91164 

8317.2 

ti  15 

12S1.8 

ti  65 

4999.4 

tms 

4389.2 

1 1 165 

83SS.1 

ti  16 

1299.9 

ti  66 

3973.6 

till! 

4439.2 

81166 

8399.9 

ti  17 

1297.1 

ti  67 

3951.6 

till? 

4471.1 

8H67 

■446.9 

ti  11 

1426.1 

tl  61 

4991.3 

tma 

4512.9 

ti  168 

I4M.9 

ti  19 

1426.2 

tl  69 

4961.3 

tni9 

4553.9 

tH69 

8521.1 

tl  M 

1S6S.9 

ti  79 

4291.9 

ti  129 

5465.3 

H 179 

•562.1 

t*  21 

1559.S 

li  71 

4184.2 

tn2i 

5596.3 

H171 

■693.7 

tl  22 

1699.2 

ti  72 

4323.1 

tn22 

5547.2 

H172 

•144.7 

ti  23 

1617.7 

ti  73 

4299.2 

t«123 

5518.2 

H173 

•615.6 

tl  24 

1827.3 

li  74 

4431.9 

tl  124 

5629.1 

H 174 

■726.6 

tl  25 

1812.3 

ti  76 

4413.5 

tn2S 

6679.9 

hits 

•767.5 

tl  26 

1952.9 

ti  76 

4553.2 

IH26 

6711.9 

HITS 

MM.  4 

ti  27 

1933.7 

li  77 

4527.2 

ti  127 

6751.9 

8U77 

N49.4 

ti  21 

2973.3 

tl  71 

4666.1 

ti  128 

6792.1 

tH78 

•899.3 

ti  29 

2952.2 

ti  79 

9196.5 

tH29 

5133.1 

tn79 

•931.2 

tl  39 

2931.9 

ti  99 

•611.4 

IH39 

5174.7 

HIM 

•972.2 

ti  31 

2171.6 

t>  11 

mi.  5 

ti  131 

5915.7 

tun 

9913.1 

ti  31 

•162.2 

ti  K 

1474.4 

tH32 

5956.6 

tni2 

9954.1 

ti  33 

2391.1 

tl  13 

•649.1 

tn33 

5997.5 

tun 

9996.9 

9i  34 

2212.1 

tl  14 

1391.4 

8H34 

6931.5 

H1I4 

9135.9 

ti  36 

2422.6 

tl  15 

•474.9 

IM35 

7129.9 

mis 

9176.9 

tl  36 

2491.7 

ti  16 

1157.4 

tU36 

7179.9 

tH86 

9217.1 

ti  37 

2641.3 

ti  17 

•321.2 

tl  137 

7211.1 

tni7 

9251.7 

tl  3S 

2611.1 

ti  m 

■939.7 

tU31 

7252.7 

till! 

9299.7 

tl  39 

2667.1 

ti  19 

1291.7 

tH39 

7293.7 

tnw 

9349.1 

tl  49 

2632.9 

ti  99 

7943.2 

ti  149 

7334.6 

81199 

9381.6 

tl  41 

2772.5 

It  91 

9119.9 

H141 

7375.6 

ti  191 

9422.5 

tl  42 

2746.2 

ti  92 

7964.2 

tn42 

7416.5 

ttl92 

9463.4 

tl  43 

2726.S 

tl  93 

7643.5 

IH43 

74S7.4 

tH93 

9S94.4 

tl  44 

2966.2 

tl  94 

7199.3 

11144 

7491.4 

tn94 

9646.3 

tl  46 

2945.1 

ti  95 

4719.7 

tins 

7539.3 

tn95 

9616.3 

ti  46 

2914.1 

ti  96 

4797.6 

tins 

7519.3 

tn96 

9627.2 

tl  47 

2991.7 

ti  97 

4669.9 

tl  147 

7621.2 

81197 

9668.1 

tl  48 

3191.4 

tl  91 

4634.9 

tl  141 

7662.1 

tnsi 

9799.1 

tl  49 

3977.4 

tl  99 

4591.1 

II 149 

7793.1 

tn99 

9759.9 

tl  69 

3217.9 

tl  199 

4665.2 

tHS9 

7744.9 

tl299 

9799.9 

1  -  to  1 

laalyalo  a< 

tala  with  otker  acceleration  data 

I  -  lack  to  MOMS  aaln  aeavo 
(EfffE*)  to  coatlaaa 
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Table  2;  Table  of  possible  materials 


chooae  the  Material  for  the  OGIUE  part  of  the  projoctlla< 


1  Structural  Stool 

2  Carbon  Stool  (yield  atrongth  •  33600  pal) 

3  Silicon  Stool  (yiold  atrongth  *  45000  pal) 

4  Nickel  Stool  (yiold  atrongth  ■  SSH4  pal) 

5  High-Strength  Steel 

6  Low  Carbon  and  Low  Alloy  Stool 

7  Coat  Iron 

8  Structural  Aluninuo  606 l -T6  or  6062-T6 

9  Structural  Atuoinuo  2014-T4 
1#  Structural  Aluolnuo  2024-T3 

11  Structural  Aluolnuo  2024-T4 

12  Structured  Aluolnuo  7075-T6 

13  Structured  Nagnoaluo  Alloy  AHC  585-T51 

14  Structured  Nagnoaluo  Alloy  AHC  58S 

15  Structured  Nagnoaluo  Alloy  AHC  S75 

16  Structured  Nagnoaluo  Alloy  AHC  525 

17  Structured  Nagnoaluo  Alloy  AH  3S 

18  Other  eater la  la 

99  lack  to  HAIN  NENUE>  2 
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APPENDIX 


EQUATIONS  FOR  MAXIMUM 
ALLOWABLE  DYNAMIC  STRESS 


L 

Material  R 


Structural  L  . 

Steel  R  <  Cc 

C  <  £  <  200 
c  R 


£  •  ^  *  allowable  unit  load  -7^-2 
A  in 


Q  149,000,000 

1  a)2 


Al 


1 


Nickel  Steel 


High-Strength 

Steel 


APPENDIX  (Continued) 


t>12° 


£  _ _ 30,000 _ 

A'  u.a!K 


0  <  ~  <  140 


140  <  £  <  200 


0  <■  |  120 


120  <  £  <  200 


0  <  £  <  110 
ft 


110  <  |  <  200 


f-  15.000  -  .325  (i) 


4. 

A 


15,000 


.5  + 


— L-.  (k 

15,860  k  R 


|  -  20,500  -  .605  (j'j 


A 


20,500 


5  + _ I _  (k\ 

°  11,630  VR/ 


2  /lY 

A  -  22,500  -  .738  (JJ 


2  .  22,500 

A  " 


.5  + 


— J _  /iY 

10,460  V  R/ 


for  y  -  33K 


for  y  -  33K 


y  “  45K 


y  -  45K 


y  -  50K 


y  -  50K 


A2 


APPENDIX  (Continued) 


0  <  j  <  105 


3-  -  25,000 
A 


^  y  -  55K 


105  <  |  <  200 


a 

A 


.5  + 


25,000 

1 

9,510 


^y  -  55K 


Low  Carbon  &  ^ 

Low  Alloy  Steel  ^  1®! 


^  4  135 
R 


2  -  36.000  -  1.172  (5^55)  forS- 
2. 79.500  -  51. 9(^35)  f°lV7 


£  <  110 
R 


a 

A 


113,000  -  11.15 


for  ^y 


h  <  95  ^  ■  145,000  -  18.36^  for^y 

i  <  90  5  *  179,000  -  27.95  (5^5) 

& 


A3 


36K 


-  75K 


-  103K 


-  132K 


-  163K 


APPENDIX  (Continued) 


Cast  Iron 


i<i°° 


£  -  12,000  -  60  \ 
A  R 


i<7° 


2 

A 


9,000  -  AO 


L 

R 


Structural 

Aluminum 

6061- T6 

6062- T6 


—  <  10 
R 


£  -  19,000 
A 


10  <  |  <  67 


£  -  20.A00 
A 


L 

R 


>  67 


£ 

A 


51.000,000 

a)2 


Structural 

Aluminum 

2014-T4 


£ 

A 


F 

co 


i.imr 


.  j  1.5E 

'  /  F 
y  co 


where  7 


l1  +  700.000  ^ 


APPENDIX  (Continued) 


and 


Structured 

Aluminum 


Structured 
Magnesium  Alloy 


F  -  35,000  for  2014-T4 
cy 

F  -  42,000  for  2024-T3 
cy 

F  -  40,000  for  2024-T4 
cy 

F  -  35,000  for  606I-T6 
cy 


^  <  1.414Y  /  ^  -  F 

R  V  Fco  A  C° 


where  F  “  1.075  F  and  F  “  66,000  for  7075—T6 
co  cy  cy 


£  - . 

to  exceed 

A 

W? 

not 
l2 
'  R2 

ALLOT 

jL 

.si 

AMC585-T51 

160,900 

.00249 

36,000 

AMC585 

46,000 

.00072 

22,000 

AMC575 

34,300 

.00053 

19,000 

AMC525 

25,500 

.00040 

16,000 

AM35 

16,750 

.00026 

11,000 
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APPENDIX  (Concluded) 


For  ocher  material  use  the  following: 


a*  L  ^  30  Then  the  max.  allowable  stress  Is  equal  to  the  yield- 
R  point  stress  of  material,^  y. 


b.  30<  L  /  100 
R  S‘ 

(^T  )  static  - 
co 


Then  max.  allowable  stress  is  given  by: 

/  > 

/fei_  fr7T~  \ 


1  +  .25  sec 


where  Cc0  «  critical  buckling  load,  lb. 

1  -  least  moment  of  inertia  of  cross  sectional  area,  in* 

2 

A  »  cross  sectional  area,  in 

R  -  least  radius  of  gyration  of  cross  sectional  area(B« 


b' 


yield-point  stress  of  material,  PSI 


L  “  length  of  column,  in 
E  -  modulus  of  elasticity,  PSI 


«  i>  100  then  (5-  >  -  $h_ 

* '  °  M2 


A6 


>|m 


